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(57) Abstract: The invention 

concerns an electromagnetic 
simulation algorithm, which enables 
to calculate the electromagnetic wave 
diffracted by a conductor operating 
in single-frequency. The invention 
concerns an electromagnetic simulation 
algorithm based on iterative resolution 
of a system of integral equations 
comprising a preconditioner Said 
preconditioner is in particular derived 
from the adaptation of Calderon 
formulae to integral equations of the 
electromagnetic boundary, also known 
as Electric Field Integral Equation 
(EFEE). It further consists in using an 
original representation of the residue of 
the calculations at each iteration. Said 
representation, as well as a projection 
and a composition are involved in the 
expression of said precondition factor. 
The invention is in particular applicable 
to simulation instruments used for 
designing receiving and transmitting 
antennae such as antennae of cellular 



ON telephones, of anti -collision radar, of 

electronic countermeasure system, of surveillance or tracking radar, or of satellite. The invention is also useful for calculating radar 
1^ reflective area of objects whereof the geometrical properties are known. 
ON 

(57) Abreg£ : La presente invention se rapporte a un algorithme de simulation 61ectromagn£tique, qui permet de calculer l'onde 
O* 61ectromagn6tique diffractee par un conducteur en regime monofrequentiel. A cet effet, Finvention concerne un algorithme de si- 

Omulation 61ectromagnetique bas6 une resolution iterative d'un systeme d'equations int6grales comprenant un preconditionneur. Ce 
preconditionneurest notamment issu de l'adaptation des formules de Calderon aux equations integrates 
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de frontiere de F61ectromagn6tisme, encore connues sous le nom de Electric Field Integral Equation (EKE) dans la literature anglo- 
saxonne. On utilise aussi une representation origniale du residu des calculs lors de chaque iteration. Cette representation, ainsi 
qu'une projection et une composition interviennent dans l'expression dudit preconditionneur. Elle s'applique notamment aux outils 
de simulation utilises lors de la conception d'antennes de reception ou d' emission telle que les antennes de telephone cellulaire, de 
radar anti-collision, de systeme de contre mesure ^lectronique (CME), de radar de veille ou de poursuite, ou de satellite. L' invention 
s'applique aussi au calcul des surfaces equivalentes radar (SER) d'objets dont on connait les propri6t£s g^ometriques. 
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ALGORITHME DE SIMULATION ELECTRO M AGNETIQUE , 
NOTAMMENT DES PERFORMANCES D'UNE ANTENNE. 

La presente invention se rapporte a un aigorithme de simulation 
electromagnetique, notamment des performances d'une antenne, qui permet 
de calculer I'onde electomagnetique diffractee par un conducteur en regime 
monofrequentiel. Elle s'applique notamment aux outils de simulation utilises 
5 lors de la conception d'antennes de reception ou d'emission telle que les 
antennes de telephone cellulaire, de radar anti-collision, de systeme de 
contre mesure electronique (CME), de radar de veille ou de poursuite, ou de 
satellite. L'invention s'applique aussi au calcul des surfaces equivalentes 
radar (SER) d'objets dont on connaTt les proprietes geometriques. 

10 On utilise les simulations d'antennes pour limiter le nombre de 

maquettes et de prototypes lors de la conception desdites antennes. Ces 
simulations permettent notamment de calculer le diagramme de 
rayonnement en champ lointain des antennes et d'adapter les antennes en 
emission ou en reception, en presence ou non d'une structure environnante. 

15 Eiles utilisent comme donnee d'entree un maillage de Tantenne dont on 
souhaite evaluer les performances, ainsi que les caracteristiques de 
Texcitation electromagnetique a laquelle elle est soumise. L'invention ne se 
limite pas aux simulations d'antennes. Elle s'applique aussi par exemple aux 
calculs de SER de cibles. On decrit maintenant, a titre d'illustration, 

20 I'application de Tinvention dans les simulations d'antennes, en reception. 

On distingue deux principles methodes dans les simulations 
couramment employees. Une premiere m6thode est basee sur le calcul par 
differences finies, encore connue sous le nom de methode d'elements finis 
. de volume. Selon cette methode, on utilise un maillage d'un volume 

25 entourant Tantenne. Un inconvenient de cette methode est que le maillage 
est necessairement borne, alors cju'on s'int§resse au diagramme de 
rayonnement a Tinfini. On doit alors realiser un compromis entre la dimension 
du volume maille, c'est a dire la precision de calcul, et le temps de calcul. 
Pour pallier cet inconvenient, on utilise une seconde methode basee sur des 

30 equations integrates dans le domaine fr6quentiel. Selon cette m6thode, on 
utilise un maillage en surface de I'antenne uniquement. On calcule 
directement le diagramme de rayonnement a I'infini a partir de courants 
electriques et magn6tiques sur la surface de Tantenne. 

Certaines techniques connues utilisant les equations integrates 

35 calculent les courants electriques et magnetiques (desquels on deduit le 
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diagramme de champ rayonne a I'infini) grace a une factorisation d'une 
matrice. Cette matrice est connue sous le nom de matrice d'interaction ou 
encore matrice d'impedance. Cette factorisation permet un calcul direct, c'est 
a dire non iteratif des courants de surface. Un inconvenient de ces 

5 techniques est que le temps de calcul est long. Si on note N le nombre de 
points intervenant dans le maillage de I'antenne (encore appele triangulation 
de surface), le temps de calcul selon ces techniques varie comme N 3 . Or le 
nombre de points N est lui-m§me lie a la longueur d'onde (et par consequent 
la frequence) de I'onde rayonnee par I'antenne. Supposons que Ton realise 

10 une simulation a 10GHz, en utilisant N points dans le maillage de Tantenne, 
et que le temps de calcul est x. Pour transposer cette simulation a 20GHz, il 
sera necessaire d'utiliser un maillage comportant 4xN points, ce qui 
representera un temps de calcul de I'ordre de 4 3 xx. Un probteme de temps 
de calcul se pose aussi lorsqu'on cherche a simuler des geometries 

15 complexes d'antennes, telles que les r6seaux de petites tailles. Ceci rend 
ces techniques inutilisables notamment dans les outils de conception qui 
requierent un temps de calcul reduit pour permettre aux concepteurs de 
realiser plusieurs tests. 

D'autres techniques connues utilisant les equations integrates 

20 permettent de diminuer le temps de calcul grace a une methode de 
resolution iterative. Si on note IT le nombre d'iterations, le temps de calcul 
selon ces techniques varie comme ITxN 2 . Un probleme de ces techniques 
est que rien ne garantit la convergence des calculs. En d'autres terrnes, il 
existe des formes d'antennes pour lesquelles on ne peut pas calculer le 

25 diagramme de champ rayonne avec ces techniques. 

Un but de I'invention est de pallier les inconvenients precites, et 
notamment de reduire les temps de calcul. 

A cet effet, invention concerne un algorithme de simulation des 
performances d'une antenne bas6 sur une resolution iterative d'un systeme 

30 d'equations integrates comprenant un preconditionneur. Ce preconditionneur 
est notamment issu de I'adaptation des formules de Calderon aux equations 
integrates de frontiere de I'^lectromagnetisme. En particulier dans le cas 
d'une antenne entierement metallique, on propose un preconditionneur pour 
liquation appelee « Electric Field Integral Equation » (EFIE) dans la 

35 litterature anglo-saxonne. On utilise aussi une representation originate du 
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residu des calculs lors de chaque iteration. Cette representation, ainsi qu'une 
projection et une composition interviennent dans I'expression dudit 
preconditionneur. 

L'invention a les principaux avantages suivants : 
5 ■ elle converge rapidement ; 

■ elle permet de simuler des geometries arbitrages et des excitations 
quelconques ; 

- son conditionnement est independant de la finesse du maillage ; 

■ elle s'adapte aux algorithmes bases sur le calcul d'une matrice 
10 d'impedance, en reutilisant ladite matrice d'impedance ; 

■ elle permet de traiter des antennes comportant, en plus du metal, des 
materiaux dielectriques. 

D'autres caracteristiques et avantages de l'invention apparattront 
plus clairement dans la description qui va suivre et dans les figures annexees 
15 qui represented : 

- la figure 1 , un maillage d'une antenne ; 

- la figure 2, une vue en coupe du maillage de la figure 1 ; 

- la figure 3, un detail du maillage de la figure 1 sur lequel est 
represents un champ de vecteurs ; 

20 - la figure 4, un diagramme fonctionnel d'un algorithme iteratif ; 

- la figure 5, une fonction de base constante par triangles sur un 

maillage ; 

- la figure 6, une fonction de base affine par triangles et continue 
sur un maillage ; 

25 - les figures 7 et 8, deux illustrations des performances de 

I'algorithme selon l'invention compare aux techniques connues. 

On se refere maintenant aux figures 1 et 2 qui represented un 
exemple de forme d'antenne dont on cherche a determiner le champ diffracte 
30 lorsque I'antenne est eclairee par une onde incidente. En d'autres termes, on 
cherche a simuler une antenne en reception. L'antenne prise dans cet 
exemple est une cavite spherique de demi-angle d'ouverture n/4. Le rayon 
interieur est 7/8, le rayon exterieur est 9/8 (unite de longueur arbitraire). La 
surface de l'antenne est notee r. Cette surface r est maillee par des 
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triangles. La surface r est supposee etre celle d'un conducteur parfait 
(rantenne) Q- plong§ dans le vide £2+. 

Dans cet exemple, Tonde electromagnetique incidente qui illumine 
Tantenne est une onde plane monofrequentielle. Cette onde 
5 electromagnetique incidente, de nombre d'onde k inc connu, est representee 
par deux champs de vecteurs notes E' nc et H inc correspondant 
respectivement aux champ electrique et au champ magnetique. Bien 
entendu, Tinvention ne s'applique pas aux seulement aux ondes planes. II est 
possible de substituer cette onde incidente par le champ emis par un dipole 

10 rayonnant par exemple (mode emetteur de Tantenne). 

On cherche a determiner Tonde electromagnetique diffractee par 
Tantenne a Tinfini, c'est a dire le diagramme de rayonnement en champ 
lointain. Cette onde electromagnetique diffractee est representee par deux 
champs de vecteurs notes E d,f et H d,f correspondant respectivement aux 

1 5 champ electrique et au champ magnetique. 

A partir du champ de courant ctrculant a la surface de ladite 
antenne, encore appele courants electriques et magnetiques de surface, il 
est possible de calculer le champ electromagnetique rayonne en tout point 
de Tespace. L'expression asymptotique a Tinfini du champ electromagnetique 

20 rayonne est le diagramme de rayonnement que Ton cherche a determiner. 
Ce calcul bien connu en electromagnetisme est rappele dans le document 
« Integral Equation Methods In Scattering Theory » de D. Colton et R. Kress 
- John Wiley & Sons, New-York, 1983. 

Le champ electromagnetique satisfait les equations de Maxwell 

25 dans le vide Q+ qui s'ecrivent : 

rot(E)= ioojLiH (1) 
rot(H)= -\cdzE (2) 



30 



35 



Les termes utilises dans ces Equations represented : 
rot Toperateur rotationnel ; 

E = E inc + E dlf le champ electrique total en notation complexe ; 
H = H ,nc + H dif le champ magnetique total en notation complexe ; 
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■ <d la pulsation de I'onde electromagnetique ; 

■ ii la permeabilite magnetique ; 

■ s la permittivity 6lectrique. 

On rappelle que le nombre d'onde est relie a la pulsation 
5 simplement par la relation suivante : 



(3) 



Le champ electrique de I'onde electromagnetique diffractee 
1 o E df s'exprime en tout point Bde Q+a partir des courants de surface notes u 
par les relations suivantes : 



E*(B)=k jG k (A,B)u(A)dS A +igrad B ( jG k (A,B)div A (u)dS i 

A<=r K Wr 



J 



(4) 



1 0 ikAB 

kV 1 4n AB 



(5) 



20 



25 



30 



Les termes des relations (4) et (5) represented : 
div A (u) la divergence prise en un point A du champ de vecteur Q ; 

grad B le gradient pris en un point B ; 

dS A un element differentiel de surface ; 
G k la fonction de Green standard; 
AB la distance entre les points A et B ; 
k la norme du nombre d'onde defini par la relation (3). 

L'homme de I'art pourra aller chercher d'autres elements 
techniques relatifs a ce calcul dans le document susmentionne dans la 
mesure ou celui-ci fait partie integrante de la description. 



Le calcul des courants de surface, c'est a dire du champ de 
vecteurs u , est determine a partir de I'equation variationnelle suivante : 



Vv m(u,v) = l(v) 



(6) 
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dans laquelle 

m(u,v)=k fjG k (A f B)u(A) v(B)dS A 6S t 



AeT 
BeT 



- 1 JjG k (A,B) div A (u) div B (v) dS A dS B (7) 

K Aer 



BeT 



et 

l(v)=- |E inc (A).v(A)dS A (8) 

Aer 

10 Cette equation variationnelle (6) est connue sous le nom de 

d'equation integrate de frontiere de relectromagnetisme, ou encore « Electric 
Field Integral Equation » (EFIE) dans la litterature anglo-saxonne. 

Afin de r6soudre cette Equation variationnelle (6) dans une 
15 simulation numerique d'antenne, on doit approcher la solution u dans un 
espace de dimension finie, dit espace de discretisation. Cet espace contient 
des champs de vecteurs qui represented des courants de surface. La 
dimension cet espace est le nombre de composantes necessaires pour 
decrire entierement ledit champ de vecteurs, en tout point de la surface r. 
20 Cette surface etant maillee, le nombre de composantes servant a decrire Ci 
sera fonction notamment du nombre de points du maillage N, ainsi que de la 
nature des fonctions de base servant d decrire le champ de vecteurs (par 
exemple des fonctions lineaires ou de degre 2). Dans la suite de la 
description, on prendra a titre d'illustration Tespace de Raviart-Thomas de 
25 plus bas degre sur le maillage de r. Cet espace de discretisation, qui 
represente des courants de surface, est note Dh. Le symbole h represente la 
longueur caracteristique, ou encore precision, du maillage. En efTet, le 
nombre de dimension de cet espace Dh depend du nombre de point N du 
maillage, qui lui-meme depend de la precision h du maillage. 
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On se refere maintenant a la figure 3 pour decrire une base de 
I'espace Dh. Dans cette base, le champ de courants de surface u est 
represents par les coefficients d'un vecteur note U. 

La base de I'espace des courants de surface contient des 
5 elements notes <Pj ou i est un indice entier associe a une arete du maillage 

de la surface r. Ces elements sont des champs de vecteurs definis sur le 
maillage de la surface r. Le champ de vecteurs <p, t represents par des 

Heches sur la figure 3, a un support born§ a deux triangles Ti et T 2 du 
maillage. Ces triangles Ti et T 2 ont en commun I'arete d'indice i et de 
10 longueur i x . Cette arrete est orientee par un vecteur de norme unitaire note 

y, . On note Pi le sommet du triangle Ti non contenu sur I'arete i ; on note P 2 

le sommet du triangle T 2 non contenu sur I'arete L On note Si et S 2 les 
surfaces des triangles Ti et T 2 . Soient Z 1 et z 2 les vecteurs de norme 

unitaire, ayant une direction normale a la surface des triangles Ti et T 2 , et 
15 orientes de I'interieur Q- vers Texterieur Q+. On defini les vecteurs x!, et 

x^les vecteurs de normes unitaires tels que le triplet (x^y^z.,) et le triplet 

(x2;y,;z 2 ) soient des triedres directs. On dSfinit le champ de vecteurs <p, 

pour tout point A appartenant a la surface r avec les relations suivantes : 

20 ■ si A e 9|(A)= ±— L P 1 A avec 9,(A)*x\>0; (9) 

1 

■ si A € T 2 , 9,(A) = ±~^A avec 9j(A)- > 0 ; (10) 

S 2 

■ sinon, 9,(A) = 6. (11) 

Une telle base de I'espace des courants est connue sous le nom 
25 de base usuelle de I'espace Raviart-Thomas ou encore courants 
elementaires de Rao-Wilton-Glisson. L'homme de I'art pourra aller chercher 
d'autres elements techniques dans « Electromagnetic scattering by surfaces 
of arbitrary shape » de S.S.M. Rao, D.R. Wilton et A.W. Glisson - IEEE 
Trans. Ant. Prop. AP-30, pp. 409-418, 1982 - dans la mesure ou ce 
30 document fait partie integrante de la description. 
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On note U* les coefficients du vecteur U. Le champ de vecteurs 
u(A) utilise dans la relation (4) se decompose sur la base usuelle de 

i'espace de Raviart-Thomas de la maniere suivante : 
5 u(A) = 2U i 9 i (A) (12) 

i 

On transpose maintenant en notations matricielles les relations 
(6), (7) et (8) en utilisant la base ((ft) decrite ci-dessus. L'equation 

variationnelle (6) s'ecrit sous la forme du systeme d'equations lineaires 
10 suivant: 

M U = L (13) 

Les termes de la relation (13) sont les suivants : 
15 ■ M est une matrice connue, dite matrice d'interaction ou encore matrice 
d'impedance ; 

■ L est un vecteur connu, dont les coefficients represented I'onde 
incidente, c'est a dire I'excitation electromagnetique ; 

■ U est le vecteur que Ton cherche a determiner, dont les coefficients 
20 represented les courants de surface. 

On definit les coefficients Mjj de la matrice d'interaction M et les 
coefficients U du vecteur L representant I'onde incidente par les relations 
suivantes : 

25 

M, =k JjG k (A,B)(9 i (A) cp j (B))dS A dS B 

AeT 

B<=r 

-Z fjG k (A,B)div A (9,)div B (9 j )dS A dS B (14) 

K AeT 
Ber 



30 



L, = j(E^(A).cp,(A))dS A 

AeT 



(15) 
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L'homme de I'art pourra aller chercher d'autres elements 
techniques relatifs au calcul de ces coefficients dans « Approximation par 
elements finis de surface de problemes de diffraction des ondes 
electromagnetiques » de A. Bendali - These de I'Universite Paris VI, 1984 - 
5 dans la mesure oCi ce document fait partie integrante de la description. 

On se refere maintenant a la figure 4 ou est illustre un algorithme 
iteratif de resolution du systeme d'equations lineaires (13) base sur la 

10 technique de gradient conjugue. II est a noter que la technique de 
preconditionnement que nous illustrons dans le cas de I'algorithme du 
gradient conjugue s'applique aussi bien d'autres algorithmes iteratifs. On 
peut citer notamment les algorithmes du Residu Minimum Generalise, 
Gradient Bi-Conjugue, Residu Quasi-Minimum et du Gradient Bi-Conjugue 

15 Stabilise encore connus dans la litterature anglo-saxonne sous les noms de 
« Generalized Minimum Residual » (GMRES), « BiConjuguate Gradient » 
(BiCG), « Quasi-Minimal Residual » (QMR), et « BiConjugate Gradient 
Stabilized » (Bi-CGSTAB). L'homme de I'art pourra aller chercher d'autres 
elements techniques sur les methodes iteratives dans « Templates for the 

20 Solution of Linear Systems : Building Blocks for Iterative Methods » de R. 
Barrett, M. Berry, T. F. Chan, J. Demmel, J. Donate, J. Dongarra, V. Eijkhout, 
R. Pozi, C. Romine et H. Van der Vorst - SIAM (1994), Philadelphia, PA - 
dans la mesure ou ce document fait partie integrante de la description. 

La figure 4 illustre un algorithme 40 qui prend en entree la matrice 

25 M et le vecteur L de I'equation (13) et donne en sortie le vecteur U. L'homme 
de I'art pourra aller chercher d'autres elements techniques relatifs a la 
resolution de systemes d'equations lineaires dans « Iterative Methods for 
Linear and Nonlinear Equations » de C.T. Kelley - SIAM Frontiers in Applied 
Mathematics, Philadelphia, 1995 - dans la mesure ou ce document fait partie 

30 integrante de la description. 

Une premiere etape d'initialisation 41 permet d'initialiser quatre 
suites de vecteurs notes U[n], R[n], S[n] et P[n] ou n est un entier. La 
premiere suite U[n] est une solution approchee qui converge vers la solution 
recherche U. La seconde suite R[n], connue sous le nom de residu, 

35 converge vers le vecteur nul. La troisieme suite S[n], que Ton nomme residu 
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10 



25 



preconditions, converge aussi vers le vecteur nul. La derniere suite P[n] est 
connue sous ie nom de direction de recherche. Les premieres valeurs de ces 
suites sont definies par les relations suivantes : 

U[0] = 0 (16) 

R[0]=L-MU[0] (17) 

S[0]=ZR[0] (18) 

P[0]=S[0] (19) 



La relation (16) est utilisee par defaut lorsqu'on ne connaTt pas de 
solution approchee. Une variante de cette relation consiste a prendre pour 
15 U[0] le resultat d'un calcul de courant de surface realise pour une meme 
antenne mais a une autre frequence. 

Le terme Z de la relation (18) est une matrice. Cette matrice est un 
preconditionneur pour la matrice M selon invention. On rappelle qu'un 
preconditionneur pour M est une matrice approchant I'inverse de M. Selon 
20 une variante, on peut remplacer le preconditionneur Z par I'identite. En 
d'autres termes, on peut initialiser S[0] avec R[0]. 

Une seconde etape d'iteration 42 permet de calculer les valeurs 
des suites precitees a un rang n+1 a parUr des termes de rang n. Cette etape 
d'iteration utilise les relations suivantes : 



30 



U[n + l] = U[n]+aP[n] (20) 

R[n + l] = R[n]-aMP[n] (21) 

S[n + l] = ZR[n + 1] (22) 

.S[n + 1), 



P[n + l] = S[n + l] + (Rfrl ^ 1 l 



"(RR 



-P[n] (23) 



avec 
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(MP[n],P[n]> (24) 
ou ( , ) represente le produit scalaire complexe en notation 

5 matricielle. 

Une derniere etape 43 realise un test de convergence. Ce test 
peut s'exprimer par exemple par I'inegalite suivante : 



|R 


n 




0 



4^<-n (25) 



10 

ou r\ est un seuil predetermine. En d'autres termes, on compare le 
residu preconditionne normalise par rapport au un seuil predetermine ti. 
Lorsque I'inegalite (25) est verifiee, le calcul s'acheve et on prend U = U[nj 

comme solution. Dans le cas contraire, on incremente la valeur de n et on 
1 5 retourne a I'etape 42. 

Bien entendu, il est possible d'utiliser un autre test de 
convergence pour arreter les iterations. On peut par exemple remplacer 
I'inegalite (25) par I'inegalite suivante : 

1*U<„ (26) 

L'algorithme du gradient conjugu6 dans les techniques actuelles 
est utilise sans preconditionneur, c'est a dire avec la matrice Z egale a 
I'identite. 

25 L'invention consiste a utiliser un preconditionneur basesur une 

generalisation d'une formule de Calderon. L'homme de Tart pourra aller 
chercher des 6lements techniques sur cette formule de Calderon dans 
« Mathematical Methods in Electromagnetism, Linear Theory and 
Applications » de M. Cessenat - World Scientific Publishing Co., page 89, 

30 1996 - dans la mesure ou ce document fait partie integrante de la 
description. Cette formule s'ecrit : 
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JMJ9A. + J<KJ^= % / 



(27) 



10 



Dans cette relation (27), M, J et^ sont trois operateurs sur les 

champs de vecteurs tangents, c'est a dire des applications qui a un champ 
de vecteurs tangents associent un autre champ de vecteurs tangents, et I 
est ('application identite. Jest I'operateur de produit scalaire avec la normale 
a la surface de I'antenne, et M est I'operateur associe a la matrice 
d'interaction dont on cherche un preconditionneur. Les operateurs M.Jetli. 
sont definis formellement par les relations suivantes : 



M$)(B) = W |G k (A,B)u(A)dS ; 

L Aer 



1 — ~j < 
+ £9rad B 



jG k (A,B)div A (u)dS > 

VAer 



(28) 



15 



(JQ)(B) = u(B)az(B) 



(29) 



20 



25 



30 



(<£u)(B) = Jgrad B (G k (A,B))Au(A)dS / 



Aer 



(30) 



Dans les relations (28), (29) et (30), Mu t Ju % et u sont des 
champs de vecteurs tangents. L'indice t de la relation (28) represente la 
composante tangentielle du vecteur entre crochets. B est un point de la 
surface r. z(B) est un vecteur unitaire, normal a la surface r en B oriente 

vers I'exterieur. 

La Demanderesse a constate que I'operateur J tf^J ^ etant 

compact, c'est a dire n§gligeable, I'operateur 4 J 9A J est 

approximativement un inverse de I'operateur M. Ensuite, les 

preconditionneur etant definis a une constante multiplicative pres, il est 
possible d'eliminer la constante 4. Sachant que J* = - J, le preconditionneur 
selon I'invention est defini a partir 6e J* !M J. On rappelle que 7* est 
I'operateur adjoint de J , c'est a dire I'operateur v6rifiant la relation suivante : 
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Jj*Q-V= |Q- JV (31) 

Le preconditionneur selon invention est une traduction matricielle 
de I'operateur J* M J. C'est dans cette traduction matricielle qu'apparatt 
5 Tavantage d'utiliser I'operateur J*MJ plutot que I'operateur JMJ. En 
effet, la traduction matricielle Z de I'operateur J*MJ est une matrice 
symetrique, ce qui est essentiel pour les algorithmes iteratifs, tel que le 
gradient conjugue. 

10 On decrit maintenant un exemple de preconditionneur selon 

I'invention qui permet d'accelerer la convergence de I'algorithme et aussi de 
rendre cet algorithme plus stable (on converge toujours vers la solution 
quelles que soient les conditions initiales). Ce preconditionneur est adapte 
aux problemes electromagnetiques et tire parti de la structure du probleme a 

15 resoudre. 

On definlt d'abord des espaces et leurs bases associees qui 
serviront dans la suite de la description. On a deja defini la base (cp,). 

L'espace engendre par cette base, Dh, est un espace de champs de 
vecteurs tangents. En d'autres termes, Dh est un ensemble de fonctions qui 

20 a tout point de r associent un vecteur tangent £ la surface maillee r. On 
definit aussi deux espaces de fonctions Ch et Sh. Les fonctions de ces 
espaces associent une valeur scalaire a chaque point de r. 

Ch est Tensemble des fonctions constantes par triangles, c'est a 
dire dont la valeur est constante pour tout point appartenant a un triangle 

25 donne. La base de Ch que Ton utilise par la suite est Tensemble des 
fonctions notees \|/| , qui valent 1 sur le triangle d'indice i et 0 en dehors. On 

se refere a la figure 5 sur laquelle est representee une fonction de base \|/, . 

On a represents par soucis de clart§ un maillage 50 plan et regulier. Les 
valeurs de la fonction vj/ r sont representees sur un axe 51, perpendiculaire 

30 au maillage 50. La fonction \\f { representee sur cette figure vaut 1 sur le 

triangle i et 0 en dehors. 

Sh est Tensemble des fonctions affines par triangles et continues. 
Ces fonctions ont un gradient constant sur tout triangle donne. La base de 
Sh que Ton utilise par la suite est Tensemble des fonctions notees Q { , qui 

35 valent 1 au noeud d'indice i et 0 sur les autres noeuds. On se refere a la 
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figure 6 sur laquelle est representee une fonction de base 6,. On a 

represents par soucis de clarte un maillage 60 plan et regulier. Les valeurs 
de la fonction 6, sont representees sur un axe 61, perpendiculaire au 

maillage 60. La fonction 9, representee sur cette figure est affine par 

5 triangles, et vaut 1 au noeud i. Cette fonction a un support (des valeurs non 

nulles) borne aux triangles 62, 63, 64, 65, 66, 67 qui ont un sommet 

confondu avec le noeud i. 

L'homme de Tart pourra aller chercher des elements techniques 

complementaires dans « Handbook of Numerical Analysis Vol. II, Finite 
10 Elements Methods (Part 1) » de P. G. Ciarlet - Ed. J.L. Lions, North-Holland, 

1991 - dans la mesure ou ce document fait partie integrante de la 

description. 

On definit maintenant des matrices qui serviront dans la suite de la 
description du preconditionneur. Ces matrices sont basees sur les bases 
1 5 definies ci-avant. Elles s'expriment par les relations suivantes : 



20 



25 



M1 y = j9 l (A).9 j (A)dS A (32) 

Aer 

M2 0 = /^(Ajdiv^dS* (33) 

AeT 

M3 y = J(9 j (a)az(A)) -<Pi(A)dS A (34) 



AeT 



ou z(A) est un vecteur de norme unitaire, normal a la surface r 
au point A, et oriente vers I'exterieur. 



30 



M4 a = Jvi/ j (A)0 i (A)dS A (35) 

Aer 

M5 5j = J0j(A) 0|(A) dS A (36) 

AeT 

M6, = (5(9,)) (37) 
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On decrit maintenant les etapes qui permettent de calculer le 
residu preconditionne a partir du residu dans les relations (18) et (22). Ces 
relations utilisent le preconditionneur Z que Ton decrit maintenant. Dans la 
description qui va suivre, on donne une decomposition de cette matrice Z en 

5 produit de matrices. Les matrices intervenant dans ce produit sont des 
matrices creuses ou des inverses de matrices creuses. Ainsi, selon une 
variante avantageuse de I'invehtion, on utilisera preferentiellement les 
matrices creuses plut6t que directement la matrice Z, ce qui permet 
notamment de reduire la memoire utilisee et le temps de calcul. 

10 On note maintenant R le residu, et S le residu preconditionne. R 

correspond respectivement a R[0] et a R[n+1] dans les relations (18) et (22). 
S correspond respectivement a S[0] et a S[n+1] dans les relations (18) et 
(22). 

Une premiere etape consiste a utiliser une representation mixte. 

15 Cette etape sera mieux comprise a I'aide des notations vectorielles qui 
seront ensuite transposees en notations matricielles. Le residu R correspond 
a une forme bilineaire sur I'espace Dh defini ci-avant. On note p cette forme 
bilineaire. On represente p par un champ de vecteurs r1 dans Dh de 
divergence nulle et une fonction q1 dans Ch d'integrale nulle. Le champ r1 

20 et la fonction q1 sont definis par la relation suivante pour tout champ de 
vecteurs v dans Dh : 

JM(A) • v(A) dS A + Jql(A) div A (v) dS A = p(v) (38) 

Aer AeT 

25 Le champ de vecteurs r1 ayant une divergence nulle, on peut 

ecrire pour toute fonction f dans Ch et d'integrale nulle : 

Jf(A)dlv A (R)dS A =0 (39) 

AeT 

30 Dans la suite, pour ne pas compliquer inutilement la presentation, 

on ne rappelle pas que les champs scalaires consideres sont d'integrale 
nulle. 
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30 




Les relations (38) et (39) se transposent de facon matricielle 
comme suit : 

= (40) 
0 AQ1 J (0) 

5 

Dans la relation (40), on a utilise une notation matricielle par blocs, 
dans laquelle : 

■ *M2 est la matrice transposee de M2 ; 

■ 0 est un bloc nul (vecteur ou matrice) ; 

10 ■ R1 est la repr6sentation matricielle de r1 ; 

■ Q1 est la representation matricielle de q1. 

On rappelle qu'il n'est pas necessaire d'inverser la matrice definie 
par blocs pour determiner R1 et Q1 a partir de R. En effet, cette matrice est 
une matrice creuse, c'est a dire qui contient beaucoup de tenmes nuls. La 
15 resolution d'un tel systeme est bien connue de Phomme du metier. Elle est 
rappelee dans « Handbook of Numerical Analysis Vol. II, Mixed and hybrid 
methods » pp 523-640 de J.E. Roberts et J.-M. Thomas - Ed. J.L. Lions, 
North-Holland, 1 991 . Ce document fait partie integrante de la description. 

Une seconde etape consiste a projeter r1 et q1 . Ces projections 
20 se traduisent en notations vectorielles par les relations suivantes : 

r^ = P Dh (rlAz) (41) 

q2 = p Sh (ql) (42) 



ou : 

■ p Dh est I'operateur de projection dans Dh ; 

■ Psh est I'operateur de projection dans Sh ; 

■ z est un vecteur unitaire, normal a r et oriente vers I'exterieur. 

Ces projections definies par les relations (41 ) et (42) se traduisent 
de facon matricielle par : 



WO 01/95159 



m 



17 



PCT/FR01/01715 



M1 0 YR2 > j _ f M3 R1> 
0 M5 AQ2 J " (,M4 Q1 ; 



(43) 



La resolution de la relation (43) est similaire a la resolution de la 
relation (40) dans la mesure ou les matrices M1 et M5 sont des matrices 
5 creuses. On peut noter que cette resolution est encore plus facile que celle 
de la relation (40) dans la mesure ou les matrices M1 et M5 sont de plus 
symetriques, definies positives et bien conditionnees. On determine ainsi R2 
et Q2 a partir de R1 et Q1 . 

Une troisieme etape consiste a combiner r2 et q2 en utilisant la 
10 relation suivante : 



r3 = r2 - rot(q2) 



(44) 



15 



Cette troisieme etape se traduit de fagon matricielle par 
R3 = R2-M6Q2 



(45) 



On appelle j I'application vectorielle qui resulte de la composition 
20 des trois etapes decrites ci-avant. Cette application est definie par : 



j(r) = r3 



(46) 



25 



Cette relation (46) se transpose en notations matricielles par : 

J R = R3 (47) 



oil on note J la matrice correspondant a 1'application vectorielle j, 
la matrice J etant definie par le produit suivant : 



30 



J = (1 -M6» rt 

v \0 M5J 



VMS 0 Y M1 M2\ 



0 M4 



^^2 0 J I^O, 



-1 



(48) 
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II n'est bien entendu pas necessaire de calculer directement J, car 
la decomposition (48) en produit de matrices creuses et d'inverses de 
matrices creuses est plus simple a utiliser. En d'autres termes, la matrice J 
est calculee de maniere implicite lors de son utilisation. 

5 

On note z I'application vectorielle qui correspond a la matrice du 
preconditionneur Z. Le preconditionneur z est defini a partir de I'application 
vectorielle j par : 

10 z = j* o m o j (49) 



ou : 

■ j* est I'adjoint de j ; 

■ o est I'operateur de composition d'applications ; 

15 ■ m est ('application lineaire Dh -> Dh* associee a la fomie bilineaire m. 

La relation (42) se traduit de fagon matricielle par : 

Z= l J M J (50) 

20 

On se refere maintenant aux figures 7 et 8 sur lesquelles sont 
illustrees les performances d'un algorithme avec un preconditionneur selon 
I'invention par rapport aux techniques connues n'utilisant pas de 
preconditionneur. 

25 Sur la figure 7, la courbe 70 represente la fonction 

ni s pi\ 



n h-> log 



10 



is 



n 



en Tabsence de preconditionneur. La courbe 71 



J 



represente la fonction n h-» log 



10 



(I*] 




l»s[o] 


II J 



en presence du preconditionneur 



decrit ci-dessus. 

Sur la figure 8, la courbe 80 represente la fonction 



30 nn log 



10 





n 


IP 




0 


lU 



en I'absence d'un preconditionneur. La courbe 81 
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represente la fonction n h» log 



10 





n 


»l 




0 


llJ 



en presence du preconditionneur 



decrit ci-dessus. 

On constate que si on prend un critere de convergence classique 
il = 10" 4 dans la relation (25), 50 iterations suffisent avec le 
5 preconditionneur. Avec un algorithme non pr6conditionn6, on n'atteint pas 
une precision suffisante en 200 iterations. 

Ces simulations numeriques on permis de montrer que Tutilisation 
d'un preconditionneur selon l'invention pernnet d'accelerer et de stabiliser les 
algorithmes iterates. 

10 Bien entendu, invention ne se limite pas a I'exemple utilise pour 

la d6crire. II est possible notamment d'utiliser d'autres fonctions de base ou 
un maillage different de ceux pris en exemple. 

L'invention se generalise a tout autre discretisation de I'espace 
Dh. Si on prend pour Dh un espace different de Tespace de Raviart-Thomas, 

1 5 on remplace les espaces Ch et Sh respectivement par : 

Ch = {div(u) | u g Dh}+ 1 (51) 
Sh = ^ | rot(p) g Dh} (52) 

20 

En d'autres termes, 

■ Ch est un espace d'elements finis minimal tel que la divergence des 
elements de Dh soit comprise dans Ch ; 

■ Sh est un espace d'6lements finis maximal tel que le rotationnel des 
25 Elements de Sh soit compris dans Dh. 



Une application principale de invention se trouve dans les outils 
de conception d'antenne, mais invention ne se limite pas a cette seule 
application. L'invention s'applique bien entendu aussi d tout outil de 
30 simulation bas6 sur le calcul du champ rayonne par un conducteur. On peut 
citer notamment le calcul de surfaces Squivalentes radar (SER) d'objets dont 
on connaTt les propri6tes geom6triques. 
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II est a noter aussi que la technique de precondltionnement decrite 
dans les cas d'une methode iterative s'applique aussi a d'autres methodes 
numeriques rapides. Ces methodes rapldes sont basees sur une resolution 
iterative, mais seuls les termes utiles aux produits matrice-vecteur sont 
5 calcules. Ainsi, on calcule de I'ordre de N x log(N) elements de la matrice 
d'impedance, au lieu de N 2 elements selon les techniques iteratives 
classiques. En d'autres termes, on calcule la matrice d'impedance de 
maniere implicite. L'utilisation du preconditionneur selon I'invention dans ces 
methodes rapides se fait sans difficulte. Ces methodes sont interessantes 

10 pour les objets de grandes taille, c'est a dire pour N grand. C'est le cas 
notamment pour les simulations d'antennes dites sur structure. On peut citer 
notamment les Methodes Multiples Multiniveaux et les Methodes Integrates 
Adaptatives, connues aussi dans la litterature anglo-saxonne sous les noms 
de « Fast Multilevel Multipole Methods » (FMM) et de « Adaptative Integral 

1 5 Method » (AIM). L'homme de I'art pourra alter chercher des elements 
techniques sur : 

■ les methodes rapides en general dans « Fast Solution Methods In 
Electromagnetics » de W. C. Chew, J.-M. Jin, C.-C. Lu, E. Michielssen, 
J.M. Song - IEEE Trans, on Antennas and Propagation, 45(3):533-543, 

20 Mars 1997; 

■ les Methodes Multipoles Multiniveaux dans « Multilevel Fast Multipole 
Algorithm For Electromagnetic Scattering By Large Complex Objects » de 
J. M. Song, C.-C. Lu, W.C. Chew, S.W. Lee - IEEE Trans on Antennas 
and Propagation, 45(10):1488-1493, Octobre 1997 ; 

25 ■ les Methodes Integrates Adaptatives dans « AIM : Adaptative Integral 
Method For Solving Large Scale Electromagnetic Scattering And 
Radiation Problems » de E. Bleszynski, M. Bleszynski, T. Jaroszewicz - 
Radio Science, 31(5):1225-1251, 1996 ; 
dans la mesure ou ces documents font partie integrante de la description. 
30 L'invention s'etend sans difficulte aux objets (antennes ou cibles) 

comportant des materiaux dielectriques. Dans ce cas, on cherche des 
courants electriques et magnetiques equivalents sur chaque interface. La 
matrice d'interaction entre ces courants comporte des blocs diagonaux du 
meme type que la matrice M decrite ci-dessus. Un preconditionneur pour la 
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matrice d'interaction est done obtenu en considerant la matrice diagonale par 
blocs, dont les blocs sont du type de la matrice Z decrite ci-dessus. 
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REVENDICATIONS 

1. Algorithme de simulation electromagnetique, pour determiner 
Tonde electomagnetique diffractee par un corps en regime monofrequentiel, 

5 a partir d'un maillage dudit corps et de Pexcitation electromagnetique, 
caracteris§ en ce qu'il comporte au moins : 

(a) une determination d'une matrice M, dite matrice d'interaction, dont les 
coefficients sont determines a partir du maillage dudit corps ; 

(b) une determination d'un preconditionneur Z de la matrice M, ce 
10 preconditionneur etant la traduction matricielle de Toperateur J* MJ , ou 9A 

est Toperateur associe aMj Toperateur produit vectoriel avec la normale a 
la surface dudit corps, et J* Toperateur adjoint de J\ 

(b) une determination des courants qui circulent a la surface dudit corps, par 
un algorithme iteratif de type gradient conjugue utilisant ledit 

15 preconditionneur Z, Talgorithme iteratif permettant de resoudre une equation, 
dite d'equation integrate de frontiere de Telectromagnetisme, s'ecrivant sous 
forme matricielle de la maniere suivante : 

MU = L 

20/ 

ou U est un vecteur que Ton cherche & determiner, dont les coefficients 
represented les courants de surface, et L est un vecteur connu, dont les 
coefficients represented Texcitation electromagnetique ; 

(c) une determination de Tonde diffractee par ledit corps, a partir desdits 
25 courants de surface. 

2. Algorithme de simulation electromagnetique selon la 
revendication 1 , caracterise en ce que les coefficients du vecteur U sont 
exprimes dans la base usuelle de Tespace de Raviart-Thomas. 

30 

3. Algorithme de simulation electromagnetique selon Tune des 
revendicatfons precedentes, caracterise en ce que le preconditionneur Z est 
determine de maniere implicite. 



WO 01/95159 



* 



23 



PCT/FR01/01715 



4. Algorithme de simulation electromagnetique selon I'une des 
revendications precedentes, caracterise en ce que le preconditionneur Z est 
defini par la relation suivante : 

5 Z=*J M J 

ou J est une traduction matricielle de I'operateur J, et *J la matrice 
transposee de J. 

10 5. Algorithme de simulation electromagnetique selon la 

revendication 4, caracterise en ce que la matrice J est definie par la relation 
suivante : 

, a, *>J m 0 VY M3 0 Y M1 M2VVl > | 
J = (1 MsJ U M4JU2 0 J UJ 

15 

oD 

M1y= j9i(A)-9j(A)dS A 

AeT 

20 M2 0 = J\j/ j (A)div A (9 J )dS A 

AeT 

M3 8 = j(9,(A)Az(A)).9,(A)dS A 

Aer 

M4 ij= J ¥j (A)0 i (A)dS A 

Aer 

M5,- Je j (A)6 i (A)dS A 

AeT 

M6 fJ = frtfo)^ 

30 OU 
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<Pi est la base usuelle de I'espace de Raviart-Thomas ; 



6i est la base des fonctions affines par triangles ; 
5 est la base des fonctions constantes par triangles ; 

A un point appartenant a la surface r dudit corps ; 

z(A) est un vecteur de norme unitaire, normal a la surface dudit 
1 o coprs au point A, et oriente vers I'exterieur. 

6. Algorithme de simulation electromagnetique selon Tune des 
revendications precedentes, caracterise en ce que I'algorithme iteratif utilise 
est un algorithme rapide, du type methode multipSles multiniveaux. 

15 

7. Algorithme de simulation electromagnetique selon Tune des 
revendications 1 a 5, caracterise en ce que I'algorithme iteratif utilise est un 
algorithme rapide, du type methode methode integrate adaptative. 

20 8. Algorithme de simulation electromagnetique selon Tune des 

revendications precedentes, caracterise en ce que le corps est une antenne 
dont on cherche a determiner une forme optimale, en utilisant I'algorithme de 
simulation dans un outil de conception d'antennes. 

25 9. Algorithme de simulation electromagnetique selon Tune des 

revendications 1 a 7, caracterise en ce que ie corps est un objet de forme 
connue dont on cherche a determiner la surface equivalente radar (SER). 



30 
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